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Cooper problem for interacting fermions is solved in a lattice. It is found that the binding energy 
of the Cooper problem can behave qualitatively differently from the gap parameter of the BCS 
theory and that pairs of non-zero center of mass momentum are favored in systems with unequal 
Fermi energies. 



I. INTRODUCTION 

Cooper found [l[ out that in the presence of the quies- 
cent Fermi sea attractively interacting electrons can form 
a bound state even if their total energy is larger than 
zero. In the absence of Fermi sea such bound state does 
not exist in a three-dimensional system in a continuum 
although such a bound state may appear in quasi two- 
dimensional systems [2]. The possibility a bound state 
indicates instability and and is traditionally thought of 
as an indicator of the instability towards the formation of 
a superfiuid. However, the proper theoretical description 
of the superfiuid itself involves a more refined many-body 
theory than the one used in the Cooper problem. The 
purpose of this report is to solve this classic Cooper prob- 
lem in a lattice for a polarized system. 

After formulating the Cooper problem in a lattice sev- 
eral questions are addressed. First, how does the bind- 
ing energy behave and is it always qualitatively similar 
to the gap parameter in the BCS theory 3]? Second, 
what happens to the bound state in systems with un- 
equal Fermi surfaces? Third, can bound states of non- 
zero center of mass be favorable in systems with unequal 
Fermi surfaces? This question is of interest since in the 
BCS theory one expects modulated, FFLO-type 0, @, 
order parameters in systems with mismatched Fermi sur- 
faces. Since optical lattices for ultracold atoms can be 
made anisotropic and spin-dependent, the problem is 
framed and solved in a rather general setting with spin- 
dependent lattice potentials. 

Previously some aspects of the two-body problem in a 
lattice have been investigated in the absence of Fermi seas 
in Refs. @, 0, S|. Dimensionality effects in the Cooper 
problem were discussed by Esebbag al. [9(. 



II. FORMALISM 

We assume a two-component system and we label the 
spin states by a = {T>4}- If the system is composed 
of ultracold neutral fermionic atoms, these "spin" states 
would correspond to either different atoms or different 



hyperfine levels of the same isotope. We further as- 
sume that atoms are at zero temperature and in a cu- 
bic lattice which is deep enough so that only the lowest 
band must be considered. Also, since the lattice is deep 
it is enough to consider only the leading order nearest 
neighbor tunneling processes with tunneling strengths t a 
where bold-face indicates that the tunneling strengths are 
represented as vectors and that the tunneling strength 
can be different in different directions. 

At zero temperature fermions on different spin states 
can interact via s-wave interaction and we take the inter- 
action between unequal fermions located at lattice sites 
X| and Xf to be <?<5(x^ — xj). Two-body wavefunction for 
the atoms is then a solution of the Schrodinger equation 



^(xi,x T ) = £V(xi,x T 



Here t • V 2 denotes the discrete kinetic energy operator 
and it acts as 



t a [ip(x a + d) - 2ip(x a ) + tp(x a - d)] 



where d = 1 is the lattice constant and a £ {x,y,z}. 
Since the interaction only depends on the relative coor- 
dinate it is useful to write the problem in terms of the 
wave- function for the relative coordinate x = xj_ — xj. 
The center of mass and relative motion no longer sepa- 
rate in a lattice, but interactions do not couple different 
center of mass states to each other. We write the wave- 
function as 



V>(x|,x T ) = U a&{Xt y !Z ye 



i K a ( c a x i ( a + ( 1 — c c 



)*T,«)^,( r ), (1) 



where c is yet to be determined coefficient related with 
the center of mass coordinate in a lattice. Since the center 
of mass and the relative coordinate do not separate, the 
algebra is somewhat more complicated than in a contin- 
uum, but nevertheless straightforward. It turns out that 
the coefficients c a are determined by the equation 



tj,a _ sin (K a (l - c a )) 
tf >a sm(K a c a ) 



(2) 
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If the lattice is spin independent the left hand side is 
equal to one and c a = 1/2. This corresponds to the usual 
center of mass transformation for equal mass atoms. If "f- 
component is very heavy, the left hand side becomes very 
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large and the coefficient c a — > 0, which again conforms to 
our intuition about the center of mass coordinate. How- 
ever, in case of general spin dependent lattice potentials, 
coefficients must be solved using Eq. ©. 

The Schrodinger equation for the relative coordinate 
then becomes 

E(K) Q Vi + 2 (t Ua + i T , Q + E(K) Q ) + 

+ fl <y(x)]V(x) = ^(x), (3) 
where 

E(K) Q = -t Lt a cos{K a c a ) - t ha cos(if a (l - c a )). 
By expanding the relative coordinate wave function 

^ x ) = ^E^ q ' x = ( 4 ) 

q 
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FIG. 1: Binding energy in a symmetric lattice with g = — 3t x 
as a function of Fermi energy (top) and effective mass ratio 
(bottom) with fi a = 3t a . 



where N is the number of lattice sites, we find 
1^£{[2E(K) Q (cos(q Q )-l) 

q a 

+ 2(t i , Q +t T , Q + E(K) Q )]-£}7V^ x 

= ^E^ q ' x,5 ( x )- ( 5 ) 

q 

By multiplying with e~ , summing over lattice sites, 
and defining a — ^ ipq, we find an equation 



g N^f(K,q)-E' W 

where 

/(K,q) - ^2E(K) a (cos(q a )-l) 

a 

+ 2 (t l>a + t T>a + E(K) a ) . 

The crucial ingredient of the Cooper problem is the 
presence of the occupied Fermi seas. Due to Pauli block- 
ing the two-body wavefunction cannot have amplitudes 
for those states that are already occupied. This restricts 
the sum-over the first Brillouin zone into just a sum 

over allowed states Y]^- If the Fermi energies are ep. ff and 
the free dispersions are e(q) ff = ^ a 2f ffj(l (l — cos(q a )), 
then in practice this means that one should only in- 
clude those modes q that satisfy e(q + cK) ^ > ep.i and 
e(— q+ (1 — c)K)f > e F ^ simultaneously. 

By using Eq. ([6]) one can solve for the pair energy E 
and investigate whether solutions with energies less than 
ep,[ + £f,t exist. For this reason it is helpful to write 
E = ep^+ep^— A so that A > for bound states. Due to 
the anisotropy and generally complicated structure of the 
Fermi surfaces, integrations are performed numerically. 



III. BINDING ENERGY IN SYMMETRIC 
LATTICES 

Fig. [1] shows the binding energy as a function of Fermi 
energy for a system where both components see an equal 
and symmetric lattice as well as a system where the ef- 
fective masses of the fermions are different. In the first 
case, the binding energy first rises quickly with the Fermi 
energy, but at e F = At the topology of the Fermi sur- 
face changes from closed to open. For higher values of 
the Fermi energy the binding energy decreases monoton- 
ically. Due to the particle-hole symmetry the gap pa- 
rameter of the BCS theory is symmetric with respect to 
€f = 6< (which corresponds to half-filling and maximum 
gap parameter in the BCS theory). Such symmetry is 
absent in the Cooper problem. The behavior at small 
Fermi energies, which correspond to low filling fractions, 
is similar to the usual free space Cooper problem. Fig. [1] 
also shows the binding energy for the case when fermions 
have different effective masses i.e. when ^ t^. As 
the t-component becomes heavier the binding energy is 
reduced. This behavior is qualitatively similar to the 
behavior of the BCS gap parameter and is to a good ac- 
curacy exponential in the mass ratio. 

In free space the solution to the Cooper problem pre- 
dicts a binding energy A ~ exp(— 2/A) while the BCS 
theory predicts a gap parameter Abcs ~ exp(— I/A), 
where A = n(ej?)|gr| and n(ep) is the density of states 
at the Fermi-energy. Therefore, for a similar coupling 
strength the binding energy in the Cooper problem is 
A ~ A^ cs , which is much smaller than Abcs m the 
weak coupling regime. In a lattice, A is also much smaller 
than Abcs j but generally the relation between these two 
quantities is more complicated. As the Fig. [T] demon- 
strates, the binding energy of the Cooper problem, is 
more sensitive to the structure of the Fermi-surfaces than 
the BCS gap parameter. 
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FIG. 2: Binding energy as a function of tunneling strength 
anisotropy for transitions into quasi-2D (top) and quasi-lD 
(bottom) situations. We choose ef,\, = £f,t = %t x and g = 



IV. BINDING ENERGY AND DIMENSIONAL 
CROSSOVERS 

If one makes the lattice anisotropic, one can effec- 
tively change the dimensionality of the system. In Fig. [2] 
demonstrates the effect of such changes in two cases. On 
the top part one changes the tunneling strength along 
^-direction while keeping other tunneling strengths the 
same. For small values of t z the lattice is deep in the 
z-direction and the system is effectively a set of quasi- 
two-dimensional systems. As t z /t x increases to /x/4 the 
structure of the Fermi surfaces change from open (cylin- 
der shaped) Fermi-surfaces into a closed (cigar shaped) 
ones. This change is reflected as a kink in the binding 
energy of the Cooper problem. 

The bottom part of Fig. [2] shows the binding energy 
for a system where t y and t z are varied, but kept equal 
to one another. For small values of t y and t z the system 
is effectively one-dimensional and Fermi surfaces are dis- 
connected. Then the Fermi surface is composed of two 
sheets with no intersecting points. At ty = t z = /i/8, 
the Fermi surfaces become connected, but are still open. 
The surfaces become connected so that the corners of 
the previously disconnected sheets merge first. As the 
tunneling strengths in y- and z-directions increase fur- 
ther, at t y = t z = (tx/4, the Fermi surfaces become closed 
and resemble the Fermi surfaces in a continuum in that 
sense. These changes in the Fermi surface topology are 
again reflected as kinks in the binding energy. Also in 
this case the binding energy is qualitatively different from 
the BCS gap parameter, which decreases monotonically 
as t y /t x = t z /t x increases. In the one-dimensional limit 
fluctuations are important and the simple BCS mean- 
field theory is unreliable. It is interesting that the two- 
body Cooper problem avoids some qualitative problems 
faced by the BCS theory. 



V. BOUND STATES OF NON-ZERO CENTER 
OF MASS MOMENTUM 

In a two-component fermionic gas the atom numbers 
of different components can be independently controlled 
and such strongly interacting polarized fermion gases 
have been recently studied experimentally [l^, [HI, 03 • 
in a harmonic trap. Studies of polarized fermionic gases 
have revealed the possibility of phase separation [13[ as 
well as the possibility of FFLO type orderparameters 
which break the translational symmetry [14j . Due to 
nesting these modulated order parameters are expected 
to be more prominent in lattices [l5| and in systems with 
reduced dimensionality (l6l . Il7j . 

In the BCS theory at T = the non-zero order pa- 
rameter can exist for Fermi energy differences of order 
SeF — €f,i — e_F,x ~ \/~2Abcs(8£f — 0) 18}. We have 
verified that a similar conclusion applies for the existence 
of the bound state in the Cooper problem in the lattice. 
I.e. the bound state exists if Scf is (roughly) less than 
the binding energy A at Sep = 0. 

In an unpolarized system chemical potentials are the 
same and the Fermi surfaces are also the same. Then 
the binding energy of the Cooper problem is maximized 
at zero center of mass momentum due to the reduced 
density of available low energy states for non-zero pair 
momentum. However, when the system is polarized the 
Fermi surfaces are different and this argument is not nec- 
essarily valid anymore. Then non-zero pair momentum 
might be favorable just like FFLO-type states can appear 
in the BCS theory. 

As Fig. [3] demonstrates, for systems with unequal 
Fermi surfaces the binding energy can indeed be max- 
imized at non-zero center of mass momentum. If the 
Fermi energy difference is too large the bound state does 
not exist, but for Fermi energy differences which are less 
than about A(Sep = 0) the possibility of non-zero center 
of mass pair must be taken into account. Fig.[3]also shows 
the binding energy in a quasi one-dimensional lattice. 
The binding energy of the non-zero center of mass pair is 
typically much larger in the one-dimensional system than 
a three-dimensional one. In the figure the binding ener- 
gies in three-dimensional and one-dimensional systems 
are roughly similar in magnitude only because the aver- 
age Fermi energy was lower in the quasi one-dimensional 
problem. One-dimensional system is qualitatively dif- 
ferent from the three-dimensional one in that as Fermi 
energy difference increases the existence of bound states 
persists to larger center of mass momenta. This is rem- 
iniscent of the BCS theory in one-dimensional systems 
where there is no upper critical polarization above which 
the order parameter disappears. 



VI. CONCLUSIONS 

In summary, we solved the Cooper problem in a lattice 
and found that the behavior of the binding energy of the 
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Cooper problem is qualitatively different from the behav- 
ior of the gap parameter of the BCS theory. Also, for sys- 
tems with different Fermi energies, if the instability exists 
it is towards the formation of pairs with non-zero center 
of mass momentum. In three dimensions at zero tem- 
perature the BCS theory for a polarized system predicts 
FFLO orderin g on ly above certain non-zero Fermi energy 
difference [141 Il5j . For smaller Fermi energy differences 
phase separation is expected to occur [l3|. In contrast, 
the Cooper problem with mismatched Fermi surfaces al- 
ways predicts instability towards non-zero center of mass 
pairs. 



FIG. 3: Binding energy as a function of pair center of mass 
momentum (along 2>direction) for few different Fermi energy 
differences SeF- Top figure: the three-dimensional lattice 
was symmetric, If = (cf,t + £f,i) /2 = 4t x , g = — 4t x , and 
(from top to bottom) Se F /e F = {0.04,0.05,0.06,0.07,0.08}. 
Bottom figure: lattice was quasi one-dimensional with t y = 
t z = t x /10, If = 3t x , g = —4t x , and (from top to bottom) 
5e F /i F = {0.03,0.05,0.07}. 



Two new interesting articles [19|, |20j discussing the 
two-body scattering in an one-dimensional lattice in the 
absence of Fermi seas have appeared after the submission 
of this report. 
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